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Abstract 

In this first paper, we demonstrate a theorem that establishes a first step toward proving a 
necessary topological condition for the occurrence of first or second order phase transitions: we 
prove that the topology of certain submanifolds of configuration space must necessarily change at 
the phase transition point. The theorem applies to smooth, finite-range and confining potentials V 
bounded below, describing systems confined in finite regions of space with continuously varying co- 
ordinates. The relevant configuration space submanifolds are both the level sets {£„ := V^{v )} v eR 
of the potential function Vjv and the configuration space submanifolds enclosed by the £„ defined 
by {M v := oo, w])}d6R, which are labeled by the potential energy value v, and where N 

is the number of degrees of freedom. The proof of the theorem proceeds by showing that, under 
the assumption of diffeomorphicity of the equipotential hypersurfaces {£„}„ e R, as well as of the 
{My}„ e ]R, in an arbitrary interval of values for v = v/N, the Helmoltz free energy is uniformly 
convergent in N to its thermodynamic limit, at least within the class of twice differentiable func- 
tions, in the corresponding interval of temperature. This preliminary theorem is essential to prove 
another theorem - in paper II - which makes a stronger statement about the relevance of topology 
for phase transitions. 

PACS numbers: 05.70.Fh; 05.20.-y; 02.40.-k 

Keywords: Statistical Mechanics, Phase Transitions, Topology 
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1. INTRODUCTION 



Some years ago, based on the well known fact that an Hamiltonian flow cor- 
responds to a geodesic flow on a suitably defined Riemannian manifold, a new 
explanation of the origin of Hamiltonian chaos has been proposed^, ^. With 
the aid of this "geometric viewpoint", the dynamical and geometrical signa- 
tures of phase transitions have been investigated in several models P, U, S, la B| • 
Invariably, the occurrence of a phase transition is signaled by a " cuspy" pattern 
of some curvature property of the underlying mechanical Riemannian manifold, 
whereas no particular pattern is displayed in the absence of a phase transition. 
On the basis of an heuristic argument jj], it has been conjectured that the 
observed geometric signatures of phase transitions could be the consequence 
of a change of the topology of the mechanical manifolds. After intermediate 
steps jg, J direct evidence has been given of the actual existence of topological 
signatures of phase transitions. These have been put in evidence through the 
numerical computation of the Euler characteristic (a topologic invariant) for 
the level sets {E^j-^gR of the potential function of a two-dimensional lattice 
(p 4 model through the exact analytic computation of the Euler character- 
istic of {M v = V^ 1 ((— oo, i>])}„eR submanifolds of configuration space for a 
mean-field XY mo 

These results have motivated the effort to make a leap forward by proving 
that topology changes of configuration space submanifolds (either £„ or M v ) 
are necessary for the occurrence of phase transitions, at least for a class of 
potentials of physical relevance. 

In the present paper, a result of this kind is actually proved in the form 
of a necessity theorem. However, one of its basic hypotheses is somewhat too 
restrictive - and cannot be relaxed in the present demonstration scheme - to 
directly use our Main Theorem as an evident rigorous support of our former 

n 

topological hypothesis]]^. This notwithstanding, the Main Theorem proved in 
the present paper is indispensable to prove a definitely stronger result, of a 
broad domain of applicability, given in paper II. 
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In the present paper, we prove the following theorem: 

Theorem 1. Let Vn (q h ...,q N ) : R N ^ R, be a smooth, non-singular, finite- 
range potential. Denote byY> v := (v), »6R, its level sets, or equipotential 
hypersurfaces, in configuration space. 

Then let v = v/N be the potential energy per degree of freedom. 

If for any pair of values v and v' belonging to a given interval 1^ = [vq, Vi) 
and for any N > N it is 

that is is diffeomorphic to Eats', then the sequence of the Helmoltz free 
energies {F N (f3)} N€N - where (3 = 1/T (T is the temperature) and (3 G Ip = 
(P(vo), P(vi)) - is uniformly convergent at least inC 2 (Ip) so that e C 2 (I/3) 
and neither first nor second order phase transitions can occur in the (inverse) 
temperature interval (/3(vq), (3{v\)). 

This is our first Theorem, given in Section El Now, for any given model 
described by a smooth, non-singular, finite-range potential, it is in general a 
hard task to locate all its critical points and thus to ascertain whether the 
theorem actually applies to it or not. Therefore we use Theorem 1 to prove - 
in paper II - a second theorem which, making a direct link between thermo- 
dynamic entropy and a weighed sum of the Morse indexes of the submanifolds 
M v , provides a general and stronger result about the relevance of configuration 
space topology for phase transitions. We anticipate below the formulation of 
this second theorem: 

Theorem 2. Let Vjv((?i, . . . , qri) '■ ^ N R, be a smooth, non-singular, finite- 
range potential. Denote by M v := V^ 1 ((— oo, v]), v E R, the generic subman- 
ifold of configuration space bounded by £„. Let {q^ G R JV }ie[i 1 A/'(t))] ^e the set 
of critical points of the potential, that is s.t. VV/v(g , c*' ) ) = 0, and J\f(v) be 
the number of critical points up to the potential energy value v. Let F(qi , £o) 
be pseudo-cylindrical neighborhoods of the critical points, and fii(M v ) be the 
Morse indexes of M v , then there exist real numbers A(N,i,Eo), gi and real 
smooth functions B(N,i,v,Eo) such that the following equation for the micro- 
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canonical configurational entropy S N \v) holds 



i r r N 

4 \v) = - log / d N q + V A(JV, z, e ) ft ^(M„_ 60 , 



n=l 



( details and definitions are given in Section 2 of paper II), and an unbound 
growth with N of one of the derivatives \d k S^\v)/dv k \, for k = 3,4, and 
thus the occurrence of a first or of a second order phase transition respectively, 
can be entailed only by the topological term J2iLoA(N,i,so) ft fii(M v - eo ) + 
Y^Sx +1 B(N,i(n),v-v? v \e ). 

Together, these two theorems imply that for a wide class of potentials which 
are good Morse functions, a first or a second order phase transition can only be 
the consequence of a topology change of the submanifolds M v of configuration 
space. 

The converse is not true: topology changes are necessary but not sufficient 
for the occurrence of phase transitions. As we point out in Remark the 
above mentioned works in Refs.^J and 11 , 3 provide some hints about the 



sufficiency conditions but rigorous results are not yet available. 

The reader can get a hold of the meaning of the main result of the present 
paper by reading just Section|2l Section|2]and the beginning of SectionElwhere 
a sketch of the proof of Lemma |U is given. In Section El we enunciate the Main 
Theorem, four main Lemmas (and give the short proofs of two of them), we 
give the condensed proof of the Main Theorem, we enunciate a Corollary to 
the Main Theorem and give its proof. 

Section apart from the already mentioned sketch of the proof of Lemma 

which is the core of the proof of Theorem 1, contains the most tedious and 
hard reading part of the paper which is necessary to prove the Main Theorem 
but not to understand the meaning of the Theorem itself. 

A preliminary account of Theorem 1 has been given in Ref . . 



2. BASIC DEFINITIONS 



For a physical system S of n particles confined in a bounded subset A d of 
IR d , d = 1,2,3, and interacting through a real valued potential function 
defined on (A d ) xn , with N = nd, the configurational microcanonical volume 
Q(v, N) is defined for any value v of the potential Vjv as 

n(v,N)= [ d qi ...dq N 5[V N ( qi ,...,q N ) -v) = [ — ^ , (1) 

J(A d )x" JS„ l|VVjv|| 

where <i<r is a surface element of := V^ 1 (f); in what follows Cl(v,N) 
is also called structure integral. The norm ||VVjv|| is defined as ||VVjv|| = 
Ei^i(^i^v) 2 ] 1 ^ 2 - The configurational partition function Z C (/3,N) is defined 
as 

da 



Z C ((3,N)= / dq 1 ...dq N exp[-pV N {q 1 ,...,q N )]= / dv e Pv 

J(A d ) xn JO JT.^ 

(2) 

where the real parameter /3 has the physical meaning of an inverse temperature. 
Notice that the formal Laplace transform of the structure integral in the r.h.s. 
of (J2J) stems from a co-area formula |l5( which is of very general validity (it 
holds also for Hausdorff measurable sets). 

Now we can define the configurational thermodynamic functions to be used 
in this paper. 

Definition 1. Using the notation v = v/N for the value of the potential energy 
per particle, we introduce the following functions: 

- Configurational microcanonical entropy, relative to E„. For any JVgN 
and v G K ; 

S N (v) = S N (v; V N ) = ^ logSl(Nv, N) . 

- Configurational canonical free energy. For any N G N and (3 G M, 

f N ((3) EE V N ) = 1 log Z c (/3, AT) . 
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- Configurational microcanonical entropy, relative to the volume bounded 
by H v . For any iVeN and v G R, 

S { ~\v) = S { N \v;V N ) = i log M(M;, TV) 

w/iere 

da 



M(u,iV)= / dq x ...dq N Q[V N (qi,...,q N ) - v] = / dry 

J(Ad)x« Jo is,, 1 1 V vjv | 

(3) 

with B[-] t/ie Heaviside step function; M(v,N) is the codimension-0 subset of 
configuration space enclosed by the equipotential hypersurface £„. TTie repre- 
sentation of M(v, N) given in the r.h.s. stems from the already mentioned 
co- area formula in flfjj . Moreover, S^ N \v) is related with the configurational 
canonical free energy, f N , for any N G N and v G R ; through the Legendre 
transform \l& 

-f N ((3)=M{(3-v-S { -\v)}, (4) 

V 

yielding, for any N G N and (3 G R ; 

-fN(J3)=/3-v N -Sk\v N ) (5) 
u>ztn, /or any iV G N and n G R ; 

««) = -^-(«), (6) 
and iae inverse relation, valid for any iVeN and /3 G R, 

= ■ (7) 



Finally, for a system described by a Hamiltonian function id of the kind 
H = Y^i=iPi/2 + Vzv((Zi) • • • j °a0> the Helmoltz free energy is defined by 

F N {/3; H) = -(Npy 1 log J d N p d N q ex V [-{3H(p, q)} , (8) 
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whence 

F N (P; H) = -(2P)- 1 log(7r//3) - f N ((3, V N )/(3 (9) 
with its thermodynamic limit (N — > oo and vol(A d )/N = const) 

Foo(/3) = lim F^;//) . (10) 

Definition 2 (First and second order phase transitions). We say that a 
physical system S undergoes a phase transition if there exists a thermodynamic 
function which - in the thermodynamic limit (TV — > oo andvol(A d )/N = const) 
- is only piecewise analytic. In particular, if the first-order derivative of the 
Helmoltz free energy F^iP) is discontinuous at some point f3 c , then we say 
that a first-order phase transition occurs. If the second-order derivative of the 
Helmoltz free energy F 0O (/3) is discontinuous at some point f3 c , then we say that 
a second-order phase transition occurs. 

Definition 3 (Standard potential, fluid case). We say that an N degrees 
of freedom potential Vjy is a standard potential for a fluid if it is of the form 

V N : B N C R N -> R 

n n 

v N (q) = £ *(lkW,ll) + $>A(g1) (ii) 

i^=j=l i=l 

where Bn is a compact subset ofM N , N = nd, ^ is a real valued function of 
one variable such that additivity holds, and where U\ is any smoothed potential 
barrier to confine the particles in a finite volume A, that is 



U A (q) 



ifqeA' 
+oo if q a A c , complement in M. N 
C°° function for q G A \ A' 



where A'c A and A' arbitrarily close to A C 1^, closed and bounded. Ua is a 
confining potential in a limited spatial volume with the additional property that 
given two limited d- dimensional regions of space, A 1 and A 2 , having in common 
ad— 1-dimensional boundary, U Al + U^ 2 = U\ lU \ 2 . By additivity we mean 
what follows. Consider two systems Si and S2, having N\ = n±d and N 2 = n 2 d 



degrees of freedom, occuping volumes Af and A 2 , having potential energies v\ 
and v 2 , for any (qi,...,q Nl ) e (Af) xni such that V Nl (q u ...,q Nl ) = v 1} for 
any (q Nl+1 , q Nl +N 2 ) e (A d ) xn2 such that V N2 (q Nl+1 , q Nl +N 2 ) = v 2 , for 
(qi, . . . , q 4- n 2 ) e (Af)*™ 1 x (Af) x ™ 2 let V N (q 1 , . . . , gAr 1+ Ar 2 ) = v be the potential 
energy v of the compound system S = S\ + S 2 which occupies the volume 
A d = Af U Af and contains N = Ni + N 2 degrees of freedom. If 

+ N 2 , Af U A d 2 ) = Vl (N u Af) + v 2 (N 2 , A d 2 ) + v'(N u N 2 , Af, A d 2 ) (12) 

where v' stands for the interaction energy between Si and S 2 , and ifv'/vi — > 
and v'/v 2 — > for N — > oo then Vn is additive. Moreover, at short distances \P 
must be a repulsive potential so as to prevent the concentration of an arbitrary 
number of particles within small, finite volumes of any given size. 

Definition 4 (Standard potential, lattice case). We say that an N degrees 
of freedom potential is a standard potential for a lattice if it is of the form 

V N : B N c R N -> R 
V N {q)= Yl C^{\\q t -q J J)+ £ $(£) (13) 

where Bn is a compact subset ofR N . Denoting by a±,...,ad the lattice spac- 
ings, if i e N d , then (iiai, . . . , iddd) £ A d . VFe denote by m the number of 
lattice sites in each spatial direction, by n = m d the total number of lattice 
sites, by D the number of degrees of freedom on each site. Thus q^ G R D for 
any i. The total number of degrees of freedom is N = m d D. Having two sys- 
tems made of N = m d D degrees of freedom, whose site indexes i^ and i^ run 
over 1 < . . . , < m, and 1 < if\ . . . , if^ < m, after gluing together the 
two systems through a common d — 1 dimensional boundary the new system 
has indexes i running over, for example, 1 < i± < 2m and 1 < i 2 , . . . , < m. 
U 

v(N + N,A d U A d 2 ) = Vl (N, Af) + v 2 (N, A d 2 ) + v'(N, N, Af, A d ) (14) 

where v' stands for the interaction energy between the two systems and if 
v' / Vi — > and v' / v 2 — > for N — > oo then Vjy is additive. 
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Definition 5 (Short-range potential). In defining a short-range potential, 
a distinction has to be made between lattice systems and fluid systems. Given 
a standard potential Vjy on a lattice, we say that it is a short-range potential if 
the coefficients are such that for any i,j elc N d , Cjj = iff \i — j\ > c, 
with c is definitively constant for N — > oo. 

Given a standard potential Vn for a fluid system, we say that it is a short- 
range potential if there exist R > and e > such that for ||q|| > Rq it is 
I^KIIllDI < IIqII where d = 1,2,3 is the spatial dimension. 

Definition 6 (Stable potential). We say that a potential Vjv is stable Jjfj/ 
if there exists B > such that 



for any N > and (q 1 ,...,q N ) G (A d ) xn , or for q} G R D , i G J C N d , 
N = m d D, for lattices. 

Definition 7 (Confining potential). With the above definitions of standard 
potentials Vn, in the fluid case the potential is said to be confining in the sense 
that it contains which constrains the particles in a finite spatial volume, 
and in the lattice case the potential Vn contains an on-site potential such that 
- at finite energy - \\qi\\ is constrained in compact set of values. 

Remark 1 (Compactness of equipotential hypersurfaces) . From the 
previous definition it follows that, for a confining potential, the equipotential 
hypersurfaces E„ are compact (because they are closed by definition and bounded 
in view of particle confinement). 

Proposition 1 (Pointwise convergence). Assume Vn is a standard, con- 
fining, short-range and stable potential. Assume also that there exists Nq G N 
such that n^>Ar dom{S { N ) and Hjv^jVo dom(SN) are nonempty sets, then the 
following pointwise limits exist almost everywhere 




(15) 



oo 



N 



lim ' (v) = Sj^ ^ (v) for v G 



N>N 




for 




N>N 
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and moreover 

oo oo 

S i - ) (v) = S 00 (v) for DG p| domiS^) n f] dom{S N ) 

N>N N>N 

Proof. The existence of the thermodynamic limit for the sequences of func- 
tions ■* and Sn, associated with a standard potential function Vn with 
short-range interactions, stable and confining is formally proved in [16j , chap- 
ters 3.3 and 3.4. To prove that in the thermodynamic limit the two entropies 
and 5*00 are equal, we proceed from the definitions of and of (3n(v), 
that is 

S ( -\v) = 1 log M(Nv,N) 

and 

noting that from the r.h.s. of Eq.(JHJ) we obtain 

dM(Nv, N) 



dv 



NQ(Nv, N) (16) 



so that 



whence 



3 (*) = 1 dM(Nv,N) = n(Nv, N) 

PN{ ' NM(Nv,N) dv M(Nv,N) 1 ' 



^\o g n(vN,N) = ^\ogM(vN,N) + ^\o g p N (v) . (18) 

Because of the existence of the thermodynamic limit j3{v) of the sequence of 
functions Pn{v) [see Proposition 2], for any given v G M. it is 

fim ^-\ogp N (v) = 

N— >oo iv 

thus, being Sn(v) = l/NlogQ(vN, N), in the thermodynamic limit, that is in 
the limit N — > oo with vol(A d ) /N = const, for any nel Eq. ()18|) implies 

S 00 (v)=S^(v) . (19) 

□ 
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Remark 2 (Equivalent definitions of entropy). In Ref.[lt] it is proved 
that the Legendre transform relating S^\v) with /jv(/3) still holds true in the 
thermodynamic limit, that is Sio\v) and /oo(^l, are related by a Legendre 



transform (see theorem 3.4-4 a ^ P-55 of Ref.\16l). Thus, after equation M9\) 



also S(v) is related with foo{P) by the same Legendre transform. 

Proposition 2 (Pointwise convergence). Assume Vjv is a standard, con- 
fining, short-range and stable potential. Assume also that there exists Nq G N 
such that HivWo dom{fjq) an< ^ VQ>n dom(j3 N ) are nonempty, then the fol- 
lowing limits exist pointwise almost everywhere 



lim f N (j3) = f(j3) , for (3 G f| dom(f N ) 

N >oo 1 1 

iV>7Vo 
oo 

lim p N (v) = P(v)) , for v G Pi dom{f3 N ) . (20) 

N — >oo 1 ' 

N>N 

Proof. See Ref.Q, chapter 3.4. □ 

Henceforth, we shall use V instead of Vn if no explicit reference the N- 
dependence of V is necessary. 



3. MAIN THEOREM 



In this Section we prove the following theorem: 

Theorem 1 (Necessity condition for Phase Transitions). Let Vn be a 

standard, smooth, confining, short-range potential bounded from below (Defi- 
nitions^^ an<20) 

V N : B N C R N -> R 
V N (g) = Yl C ii*(\\Qi-Qi\\)+ E $ ® ( 21 ) 

Let fif, $) be real valued one variable functions, let i,j label interacting pairs 
of degrees of freedom within a short-range, and let {Y, v } veR be the family of 
N — l-dimensional equipotential hypersurfaces T, v := V^ 1 (f), v G R, ofM. N . 

12 



Let Vo,V\ G M, vq < v±. If there exists N such that for any N > N and 
for any v,v' G l d = [v ,Vi\ 

Y*m is C°° — diffeomorphic to Saw, 

(notation: ^m 1 ) then the limit entropy S(v) is of differentiabil- 

ity class C 3 (Jtj) 7 and, consequently, (5{v) belongs to C 2 (Iv), whence the limit 

° ° 

Helmholtz free energy function G C (ip), where Ip denotes open interior 
of (3([v ,Vi\)), so that the system described by V has neither first nor second 

o 

order phase transitions in the inverse-temperature interval I p. 

The idea of the proof of the Theorem 1 is the following. In order to prove 
that a topology change of the equipotential hypersurfaces of configuration 
space is a necessary condition for a thermodynamic phase transition to occur, 
we shall prove the equivalent proposition that if any two hypersurfaces E^at) 
and Xv(at) with v(N),v'(N) G (v (N) , v±(N)) are diffeomorphic for all N, pos- 
sibly greater than some finite N , then no phase transition can occur in the 
(inverse) temperature interval [limjv^oo P(vo(N)), limjv->-oo fl(vi(N))]. To this 
purpose we have to show that, in the limit N — > oo and vol(A d )/N = const, 
the Helmoltz free energy F 0O (/3;if) is at least twice differentiable as a func- 
tion of (3 = 1/T in the interval [limjv^oo f3(vo(N)), limjv^oo f3(vi(N))]. For the 
standard Hamiltonian systems that we consider throughout this paper, being 
Fn(/3) = — (2/9) -1 log(ii/(3) — /W(/3)//3, this is equivalent to show that the se- 
quence of configurational free energies {/at(T; if )}Ar eN+ is uniformly convergent 
at least in C 2 so that also {/oo(T; H)} G C 2 . 

We shall give the proof of Theorem 1 through the following Lemmas, which 
are separately proven in subsequent Sections. 

Lemma 1 (Absence of critical points). Let f : M [a,b] a smooth map on 
a compact manifold M with boundary, such that its Hessian is non- degenerate. 
Suppose f(dM) = {a,b} and that for any c,d G [a, b] it is / _1 (c) ~ f~ 1 (d), 
that is all the level surfaces of f are diffeomorphic. Then f has no critical 
points, that is ||V/|| > C > 0, in [a,b]; C is a constant. 
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Proof. Since / is a good Morse function, let us consider the case of the ex- 
istence of - at least - one critical value c G [a, b] so that V/ = at some 
points of the level set / _1 (c). The set of critical points <r(c) = {x^ ki G 
/ _1 (c)|(V/)(x*' fci ) = 0} is a point set [l^], the index % labels the different crit- 
ical points and k\ is the Morse index of the 2-th critical point. After the "non- 
critical neck" theorem ^?J|, we know that the level sets f^iv) with v G [a, c—e] 
and arbitrary e > are diffeomorphic because in the absence of critical points 
in the interval [a, c — e] for any v, v' G [a, c — e], with arbitrary e > 0, f~ l {v ) is 
a deformation retraction of / -1 (t/) through the flow associated with the vector 

le neighborhood of each critical point 
allows to represent the function / 



field jig X = -V//HV/II 2 . Now, in t 
the existence of the Morse chart 



as follows 

f( X ) = f( X c % ) ~ X l ~ ' ' ' ~ X \ + X ti+1 + ' ' ' + X n ) (22) 

whence the degeneracy of the quadrics, for v = c, entailing that the level set 
/ -1 (c) no longer qualifies as a different iable manifold. Thus for any v G [a, c—e] 
and arbitrary e > 0, it is 

r\v)#r\c). (23) 

In conclusion, if for any pair of values v,v' G [a, b] one has / -1 (i>') ~ f~ l (v), 
no critical point of / can exist in the interval [a, b]. □ 

Lemma 2 (Smoothness of the structure integral). Let Vn be a stan- 
dard, short-range, stable and confining potential function bounded below. Let 
{T, v } vm be the family of (N — 1)- dimensional equipotential hyper surf aces 
'■= V^iv), v eM, ofR N , then we have: 

If for any v,v' G [uo,Ui], £« ~ £„/ then Q(v, N) G C°°(]vq, Vi[) . 



Proof. The proof of this Lemma is given in Section EJ □ 

Lemma 3 (Uniform convergence). Let U and U' be two open intervals of 
R. Let h N be a sequence of functions from U to U' , differentiate on U , and 
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let h : U — > U' be such that for any x G U, linijv^oo h^{x) = h(x). 
If there exists Mel such that for any iVeN and for any a G U it is 
A 1 < M , then h is continuous at a for any a G U . 



dx 

Proof. From the assumption that for any iV G N and for any a G U it 
is \h' N (a)\ < M, and after the fundamental theorem of calculus, the set of 



functions {/iArjjveN is equilipschitzian and thus uniformly equicontinuous 1^|. 
Then, from the Ascoli theorem on equicontinuous sets of applications 19(, 
it follows that for any a G U the closure of the set of functions {/iat}at 6 n 
is equicontinuous, and thus the limit function h is continuous at a for any 
aeU. □ 

Lemma 4 (Uniform upper bounds). Let be a standard, short-range, 
stable and confining potential function bounded below. Let {Y< v } veR be the 
family of (N — 1) -dimensional equipotential hypersurfaces T, v := V^ l {v), v G R, 
ofR N , if 



for any N, for any v,v' G h = [v ,Vi], Y, m 



then 



sup |Sjv(f)| < oo and sup 



d k S 



N 



< oo, k = 1,2,3,4. 



Proof. The proof of this Lemma is given in Section □ 

Proof (Theorem 1). Under the hypothesis that all the level surfaces of Vn 
are diffeomorphic in the interval Iy we know from Lemma H that there are no 
critical points of Vn in I v , i.e. there exists C(N) > such that for any N > Nq 

for v G Iy, and for any x G T, m , \\W N (x)\\ > C > . (24) 

Therefore, the restriction of Vn 

^ = V lv -, (Im) : V^(I m ) CB^R (25) 
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always defines a Morse function, since Vn is bounded below. Notice that 

S N (m;V N ) ]0 =S N (m;V N ) ]0 , (26) 

in what follows we shall drop the tilde and will denote the above given 
restriction. 

Now, since the condition (J24)) holds for the hypersurfaces {E^}. ° , from 

o 

Lemma El it follows that for any N > N , Q(Nv,N) is actually in C°°(/ s ), 

o 

where J B = (v ,vi); this implies that for any iV > N , also Sn belongs to 
C°°(lv). 

While at any finite N - under the main assumption of the theorem - the 
entropy functions Sn are smooth, we do not know what happens in the iV — ► oo 
limit. To know the behaviour at the limit, we have to prove the uniform 
convergence of the sequence {5jv}iveN+- Lemmas El and |U prove exactly that 
this sequence is uniformly convergent at least in the space C 3 (/ s ), so that we 

o 

can conclude that also S G C (1$). 

As S = in Iy (Proposition [TJ), also lies in C 3 (j c ) and (3 in C 2 (ly). 
Moreover, by definition and existence of the uniform limit of {«S , Ar}A r eN + , f° r 

o 

any !) 6/„- we can write 

S(v) = f(P(v))+P(v)-v 

which entails / G C 2 ((3(°h)) = C 2 (°Ip). 

Since the kinetic energy term of the Hamiltonian describing the system S 
gives only a smooth contribution, also the Helmoltz free energy has differ- 

o 

entiability class C (ip)- Hence we conclude that the system S does not undergo 
neither first nor second order phase transitions in the inverse-temperature in- 
terval /3 elp. □ 

Corollary 1. Under the same hypotheses of Theorem^ let {M v }„ e R be the 
family of the N -dimensional subsets M v := ((— oo, v]), «gK, ofR N . Let 
v ,Vi G R ; Vq < V\. If there exists N such that for any N > N and for any 
v,v' G h = [v ,vx\ 

M m is C°° — diffeomorphic to M m > , 

16 
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then the limit entropy >(v) is of differentiability class C 3 (Iy), and, conse- 
quently, P(v) = dS(~>/dv belongs to C 2 (Iv), whence the limit Helmholtz free 

° ° 

energy function F^ G C (Ip), where Ip denotes open interior of j3([v ,Vi])), 
so that the system described by V has neither first nor second order phase 

o 

transitions in the inverse-temperature interval I p. 

Proof. If for any v, v' G Iy = [^0,^1] it is Mm ~ Mm', then after Bott's 



"critical- neck theorem" 2l|, there are no critical points of Vn in the interval 
[wo,t>i]. As a consequence of the absence of critical points in [^0,^1], after the 
"non-critical neck theorem" [l3| for any v,v' G /« = [uo, ih.] it is £jyu ~ £jv«'. 
Now Theorem ^implies S(v) G C 3 ^), so that using Proposition^we have also 
S('(v) G C 3 (Iy). Then using equation (j3J) we have foo(P) G C 2 (J S ) and thus 

o 

Foe G C (ip), so that neither first nor second order phase transitions can occur 
in the inverse temperature interval Ip— (dS^ /dv\v = y , dS^ /dv\y = ^. □ 



4. PROOF OF LEMMA 2, SMOOTHNESS OF THE STRUCTURE 
INTEGRAL 

We make use of the following Lemma 

Lemma 5. Let U be a bounded open subset ofM, N , let ip be a Morse function 
defined on U , ip : U C M. N — > M and T = {T, v } v the family of hypersurfaces 
defined as E„ = {x G U\ip(x) = v}, then we have: 

if for any v,v' G [v ,«i], « Ej, 

then, for any g G C°°(U), / g da is C°° in ]v ,Vi[. 



Proof. To prove this Lemma we need the following Theorem |15l. 22]: 

Theorem (Federer, Laurence). Let O C MP be a bounded open set. Let 
ip G C n+1 (0) be constant on each connected component of the boundary dO 
and g G C n (0). 
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By introducing O t .t> = {x E O \ t < if)(x) < t'}, and F(v) = f^ =v ^ g da p ~ l , 
where da v ~ x represents the Lebesgue measure of dimension p — 1. 

If C > exists such that for any x G O t ,f, \ \Vif)(x)\\ > C, 
for any k s.t. < k < n, for any v G]t,t'[, one has 

d k F , 



dv k 



[ A k g da^ 1 . (27) 



with Ag = V (j&^g) ^ . 

By applying this Theorem to the function if> of the LemmaElwe have that, if 
there exists a constant C > such that for any x G O VOjVl it is || V^(:r)|| > C, 
then 

d k F f 

_ F ( t; ) = A k gda, Vw^o.tJil 

Now, under the hypothesis that for any v,v' G [i>o,ui], rs £.„/, 

we know from Lemma ^ "absence of critical points", that this hypothesis is 

equivalent to the assumption that for any v G [t> , t>i], has no critical 

points. Hence there exists a constant C > such that Vx G O vom \\Vif)(x)\\ > 

C. Furthermore, as HV^II is strictly positive, A is a continuous operator on 

d k F 

O vo v . Thus, being S„ compact, — — r is continuous on the interval ]w ,fi[, Vfc, 

air 

namely J s gda G C°°(\vq, v%[) . 

To conclude the proof of the Lemma 121 we have to use Lemma 03 taking 
if) — Viv and g = 1/||Wjv||, assuming that Vat is a Morse function and that 
|| VVjv|| is strictly positive (absence of critical points of stemming from the 
hypothesis of diffeomorphicity of Theorem 1). □ 



5. PROOF OF LEMMA H UPPER BOUNDS 

The proof of this Lemma is splitted into two parts. In part A some prelim- 
inary results to be used in part B are given, and in part B the inequalities of 
the Lemma |U are proved. 

The proof of Lemma |U is the core of the proof of Theorem 1. Thus, as the 
proof of Lemma|3]is lengthy, in order to ease its reading we premise a summary 
of it. 
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Sketch of the proof . 

In order to prove Theorem 1, we have to show that the assumption of 
diffeomorphicity among the T,m for v £ [vp, v i], entails that Soo(v) is three 
times different iable. After the Ascoli theorem 19|, this is proved by showing 
that for v £ h = [vo,Vi] and for any N, the function Sn(v) and its first four 
derivatives are uniformly bounded in N from above, that is, for any iV £ N 
and v £ [vq, V\] 



sup \Sn(v)\ < oo , sup 



d k S 



N 



< oo, k = 1,..,4. (28) 



After Definition 1 for the entropy, the first four derivatives of S^{y) are 

dqS N = (i/N)(dv/dv)n'/n , 

%s N = N[Q"/n - (n'/n) 2 } , (29) 

d 3 -s N = N 2 [n m /n - 3iw/fi 2 + 2(n'/n) 3 ] , 
d*s N = N 3 [a iv /a - An"'n'/n 2 - 3(n"/n) 2 + i2n"(n') 2 /n 3 - 6(fi'/fi) 4 ] , 

where the prime indexes stand for derivations of Q(v, N) with respect to v — 
vN. In order to verify whether the conditions ()28|) are fulfilled, we must be 
able to estimate the iV-dependence of all the addenda in these expressions for 
the derivatives of Sn- 

Being the assumption of diffeomorphicity of the equivalent to the ab- 
sence of critical points of the potential, we can use the derivation formula 
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da , . 

(30) 



where A k stands for k iterations of the operator 

W \ 1 



A{.) = V 



NV\ J WW 



A technically crucial step to prove the Theorem is to use the above formula 
(J3U|) to compute the derivatives of fl(v,N), in fact these are transformed into 
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the surface integrals of explicitly computable combinations and powers of a 
few basic ingredients, like ||VV||, dV/dqi, d 2 V/dqidqj, d 3 V/dqidqjdqt and so 
on. 

The first uniform bound in Eq.([2*%]). |Sjv(^)| < oo, is a simple consequence 
of the intensivity of Sn(v). 

To prove the boundedness of the first derivative of Sjy, we compute its 
expression by means of the first of Eqs.()29p and of Eq. (}3T)j) . which reads 



dS 



N 



dv 



1 



AV _ 2 Y, i j9 i Vd? j V&V 



iwi 



w 



da 



IWI 



(31) 



with diV = dV/dq % and i,j = 1, . . . , N, whence (with an obvious meaning of 
<■>=.) 



dS 



N 



dv 



< 



| AV 
IWII 



Eij d'Vd^V&V 



|W| 



(32) 



the r.h.s. of this inequality - in the absence of critical points of the potential 
- can be bounded from above by (see Lemma [HJ) 



Ay |) 



<iiwii 2 >e„ v^vy (liwuv 

As we have assumed that V is smooth and bounded below, and after the 
argument put forward in Remark |5J we have (| AV = (| ^liY I 

)e« < iVmaXj(| dfjV \)^ v and, as we have also assumed that V is a short range 
potential, the number of non-vanishing matrix elements dfjV is N(n p + 1) 
where n p is the number of neighbouring particles in the interaction range of 
the potential, thus (| &V&?jVd>V \) < N(n p + 1) max^d d i Vdf j Vd j V 

Moreover, the following lower bounds exist for the denominators in the 
inequality 

(l|Vy|| 2 ) E „ > N min^V) 2 )^, and (||W||V > 
iV 2 mm^Vf {d 3 Vfh v . 

Finally, putting m = maXij(\ d l VdfjVd j V C\ = m.mi((diV) 2 )-£ v and 
C2 = mnij jiidiV) 2 (djV) 2 )^, by substituting in Eq.© the upper bounds for 
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the numerators and the lower bounds for the denominators we obtain 



OS 



N 



dv 



ma Xi (j dlV Db. / 1 



n„ m „ / 1 , 



which, in the limit iV — > oo, shows that the first derivative of the entropy is 
uniformly bounded by a finite constant. This first step proves that *S'oo(^) is 
continuous. 

The three further steps, concerning boundedness of the higher order deriva- 
tives, involve similar arguments to be applied to a number of terms which is 
rapidly increasing with the order of the derivative. But many of these terms 
can be grouped in the form of the variance or higher moments of certain quan- 
tities, thus allowing the use of a powerful technical trick to compute their 
N- dependence. For example, using Eq.(j30J) in the expression for d 2 Sjy, we get 



d 2 S N 


< N 


dv 2 





(« 2 >E„ 



(a) 



+ N 



(VW-V(a)>3 



(35) 



where a = ||W|| A(1/||W||) and f_ = V/||W||. Now, it is possible to 
think of the scalar function a as if it were a random variable, so that the first 
term in the r.h.s. of Eq. (j35j) would be its second moment. Such a possibility 
is related with the general validity of the Monte Carlo method to compute 
multiple integrals. In particular, since the £„ are smooth, closed (V is non- 
singular), without critical points and representable as the union of suitable 



subsets of 



)JV-1 

*■ ? 



the standard Monte Carlo method 



24j is applicable to the 



computation of the averages (-)s„ which become sums of standard integrals 



m 



iiV-l 



. This means that a random walk can be constructively defined on 



any £„, which conveniently samples the desired measure on the surface (see 
Lemma 6). Along such a random walk, usually called Monte Carlo Markov 
Chain (MCMC), a and its powers behave as random variables whose "time" 
averages along the MCMC converge to the surface averages (•)■£„. Notice that 
the actual computation of these surface averages goes beyond our aim, in 
fact, we do not need the numerical values - but only the iV-dependences - 
of the upper bounds of the derivatives of the entropy. Therefore, all what 
we need is just knowing that in principle a suitable MCMC exists on each 
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£„. Now, the function a is the integrand in square brackets in Eq.()3ip. where 

the second term vanishes at large N, as is clear from Eq. (|34j) . Therefore, 

at increasingly large N, the approximate expression a = J2f=i dtiV/\\W\\ 2 

tends to become exact, a is in the form of a sum function a = iV _1 YliLi a i 

of terms a* = N8%V/\\VV\\ 2 , of 0(1) in N, which, along a MCMC, behave 

as independent random variables with probability densities itj(aj) which we 

do not need to know explicitly. Then, after a classical ergodic theorem for 

I I 

sum functions, due to Khinchin |25l |. based on the Central Limit Theorem of 
probability theory, a is a gaussian-distributed random variable; as its variance 
decreases linearly with A, limjv^oo -^|(« 2 )s„ — ( a )IJ = const < oo. 

Arguments similar to those above used for the first derivative of Sn lead 
to the result liniAr^oo N\ (ip(V) ■ ip = const < oo, which, together with 

what has been just found for the variance of a, proves the uniform boundedness 
also of the second derivative of Sn under the hypothesis of diffeomorphicity of 
the £„. 

Similarly, but with an increasingly tedious work, we can treat the third 
and fourth derivatives of the entropy. In fact, despite the large number of 
terms contained in their expressions, they again belong only to two different 
categories: those terms which can be grouped in the form of higher moments of 
the function a, and whose iV-dependence is known after the above mentioned 
theorem due to Khinchin and Lemma 7, and those terms whose iV-dependence 
can be found by means of the same kind of estimates given above for dyS^. 
Eventually, after a lenghty but rather mechanical work, also the third and 
fourth derivatives of Sn are shown to be uniformly bounded as prescribed by 
Eq.(|28|). Whence the proof of Theorem 1. 

5.1. Part A 

We begin by showing that on any (A — 1) -dimensional hypersurface T^nv — 
V^iNv) = {X e R N | V N (X) = Nv} of R N , we can define a homogeneous 
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non-periodic random Markov chain whose probability measure is the configu- 
rational microcanonical measure, namely der/||VVzv||. 

Notice that at any finite N and in the absence of critical points of the 
potential Vn (because of ||VVjv|| > C > 0) the microcanonical measure is 
smooth. The microcanonical averages ( are then equivalently computed 
as "time" averages along the previously mentioned Markov chains. 

In the following, when no ambiguity is possible, for the sake of notation we 
shall drop the suffix iV of Vn- 

Lemma 6. On each finite dimensional level set E^ = V~ 1 (Nv) of a standard, 
smooth, confining, short range potential V bounded below, and in the absence of 
critical points, there exists a random Markov chain of points {Xi G R }igN + ; 
constrained by the condition V(Xj) = Nv, which has 

as its probability measure, so that, for a smooth function F : Mr —>■ R it is 

Proof. As the level sets {Ea^I^r are compact codimension-1 hypersurfaces 
of R^, there exists on each of them a partition of unity [23]. Thus, denoting 
by {Ui}, 1 < i < m, an arbitrary finite covering of by means of domains 
of coordinates (for example by means of open balls), a set of smooth functions 
{ifi} exists, with 1 > ipi > and = 1> f° r an y point of Ejv c . Since the 

hypersurfaces E^ are compact and oriented, the partition of the unity {ifi} 
on Ejv®, subordinate to a collection {Ui} of one-to-one local parametrizations 
of Ejv s , allows to represent the integral of a given smooth (N — l)-form to as 
follows 



„ „ / m \ m „ 



^Nv \i=l / i=\ 

Now we proceed constructively by showing how a Monte Carlo Markov Chain 
(MCMC), having (J36|) as its probability measure, is constructed on a given 
Ejvtj- 
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We consider sequences of random values {xj : % G A}, with A the finite 
set of indexes of the elements of the partition of the unity on Ea^, and 



x, = (xj, . . . , x i 1 ) the local coordinates with respect to Ui of an arbitrary 



representative point of the set Ui itself. Then we define the weight tt(z) of the 
i-th element of the partition as 

-i 



7r m 



E 

,fc=i 



and the transition matrix elements 




(3? 



Pij = mm 



tt(z) 



(39) 



which satisfy the detailed balance equation n{i)pij = Tr(j)pji. Start- 
ing from an arbitrary element of the partition, labeled by io, and us- 
ing the transition probability (J39j) we obtain a random Markov chain 
{io, h ■ ■ ■ , ik, ■ ■ ■ } of indexes and, consequently, a random Markov chain of 



points {xi , Xi 17 . . . , Xi k , . . . } on the hypersurface Ejv?j- Now, let (xp, . . . 



,N-1 



, X p 



be the local coordinates of a point P on Ejv s and define a local reference frame 
as {d/dxp, . . . , d/dxp' 1 , n(P)} where n(P) is the outward unit normal vector 
at P; through the point-dependent matrix which operates the change from 
this basis to the canonical basis {ei, . . . ,eAr} of M. N we can associate to the 
Markov chain {xi , , . . . , Xi k , . . . } an equivalent chain {X io , , . . . , X ik , . . . } 
of points identified through their coordinates in M. N but still constrained to 
belong to the subset V(X) = v, that is to S^. By construction, this Monte 
Carlo Markov Chain has the probability density (j3*U|) as its invariant probabil- 



ity measure 



24j . moreover, for smooth functions F, smooth potentials V and 



in the absence of critical points, .F/||VV|| has a limited variation on each set 
U{, thus the partition of the unity can be made as fine grained as needed - 
keeping it finite - to make Lebesgue integration convergent, hence Equation 
(jHZD follows. □ 



In part B we shall need the iV-dependence of the momenta, up to the 
fourth order, of the sum of a large number of mutually independent random 
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variables. These iV-dependences are worked out in what follows by using and 
extending some results due to Khinchin j^j]. 

Definition 8. Let us consider a sequence {f]k}k=i,..,N of mutually independent 
random quantities with probability densities {uk(x)}k=i,..,N- Let us denote with 
flfc = / x Uk(x) dx the mean of the k-th quantity and with 



bk = J (x — a k ) 2 Uk{x) dx Ck = J \x — a k \ 3 u k {x) dx 
dk = (x — ak) A Uk(x) dx e-k = \x — a k \ 5 Uk(x) dx 



its higher moments. 

Theorem (Khinchin). Let us consider a sequence {rjk}k=i,..,N of mutually 
independent random quantities with probability densities {uk(x)}k=i...,N ■ With- 
out any significant loss of generality we assume that the ak are zero. Under the 
conditions of validity of the Central Limit Theorem (see E^}/ ). the probability 
density Un{x) of sn = SfcLi^fc given by 



Um(x) 



(2ttB 



exp 



N) 



X 



Sn + T N x 



R 2 

N 



+ o 



N 2 



2B N 

V I x \< 2loe 2 N 



U N (x) 



(2ttB 



exp 



N 



x 



2B 



N 



+ 



Vx e 



(40) 
(41) 
(42) 



where Bn = J2iLi h and where Sn and T/v are independent of x such that 

limjv »oo^ _1 Sn and lim^ >00 N~ x Tn are finite values (allowed to vanish) 

and where log 2 N stands for (log iV) 2 . 

Lemma 7. Consider a sequence {rjk}k=i,..,N of zero mean, mutually inde- 
pendent, random variables with probability densities {uk(x)}k=i,..,N- Denote 
with B' N , C' N and D' N the second, third and fourth moments respectively of 



'JV 



h J2k=i Vk, and with K' N = D' N — 3B' N the fourth cumulant of s' 



AT- 
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If the random quantities fulfil the hypotheses of the Central Limit Theorem, 
then 



(i) lim N B' N = est < oo 

JV — >oo 



(ii) lim N 2 C' N = 

N — >oo 

(iii) lim N 3 K' N = 

N — >oo 



Proof. Assertion (i). 

Let Bn be the second moment of sn = ^2k=i Vk- After the above reported 
Khinchin theorem, we have 



B N = 



x | 2 Un(x)oIx 
1 



(27^)3 



x I exp 



x 



2B 



N 



dx + / | x | 2 R,N(x)dx 



where Rn{x) is a remainder of order 1/JV. The r.h.s. of this equation is the 
second moment of the gaussian distribution which is just B N . Then B N can 
be rewritten, using again Khinchin theorem, as 

, 2 Sn + Tj^x 



lim Bn 

N >oo 



lim Bn + lim 

N—>oo JV^oc /| x | <2 i og 2 jV 



lim Bn + lim 

JV >oo JV >oo 



X 



X 



B 2 

JV 



= lim Bn H lim 

JV — >oo 3 N — >oo 



|x|<21og 2 JV 
JV 



i2 Sn 
B 2 



Sn \ok 6 N 



B 2 



Now let U' N (x) be the probability density of s' N = jjJ^k^Vk, its second mo- 



ment B' N is equal to 



B' 



N 



x | U N (x)dx = j-^ B N 



and thus 



lim N B' N = lim — — H lim — — r — 

JV — >oo JV — >oo N 3 N — >oo 



JV 



(43) 
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Since limjy >0 o N 1 B N is a finite non-vanishing value and lim^r >tXl N 1 Sn 

is a finite value, we conclude that 



lim N B' N = est < oo . 

N — >oo 



(44) 



□ 



Proof. Assertion (ii). 

Let Cn be the third moment of sn = XlfcLi^fc- After Khinchin theorem we 
have 



C 



N 



x | 3 Un{x)oIx 
1 



(2ttB n ) 1 2 



x | exp 



x 



IB 



N 



dx + J | x | 3 R N {x)dx 



where Rn(x) is a remainder of order 1/N. The first term of the r.h.s. is 
identically vanishing because it is an odd moment of a gaussian distribution. 
Thus Cm can be rewritten, using again Khinchin theorem, as 

,3 + Tnx 



lim Cn = lim 



N — >oo 



N »oo 



lim 

N >oo 



I 

L 



|a;|<21og 2 N 



X 



X 



5 

B~ 2 



N 



|3 



2 3 lim 

N — >oo 



S N log 8 N 



■ :2iog 2 iv B 2 N JV ^°° 5^ 

1 V^^ 



Now let U' N (x) be the probability density of s' N = jj J2 k=1 r] k , its third moment 
C' N is equal to 



C' N = I \x\ 3 U' N (x)dx =j^C N 



which leads to the conclusion 



lim tfCl,-* lim S "'° g N 



N — >oo 



= . 



N — >oo 



NBl 



(45) 
□ 



Proof. Assertion (Hi) 



Let Kn be the fourth cumulant of sn = J2k=i Vk- we have 



K 



N 



x UN(x)dx — I / x UN(x)dx 



(46) 
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which, using Khinchin theorem, can be written as 
K N — - I x a Gn{x)(Ix — ( I x 2 Gn{x)cIx 



3„ 

I r } r " 2 



+ — J x Rn(x)cIx — x R,N(x)dxj — 2 j x Rw(x)dx j x Gn(x)cIx 

1 2 

where Gjy(x) = (2ttB n )~2 exp — ^f^- is a gaussian probability distribution 
and Rn{x) the remainder of order 1/N. 

The sum of the first two terms of the r.h.s. of the equation above is the 
fourth cumulant of a gaussian distribution, thus vanishing. 

Again using Khinchin theorem we can write 

lim Kjq = - lim f x A — R N dx 

N — >oo 3 N- 



.2 

2 Sn + TffX 



lim / x g <ix 



2 S N + T N x , ., 



lim / a; g <ix / x G^(x)dx 

N ^°° J\x\<2\og 2 N B% J 

2 6 log 10 iV^ 2 8 log 12 iV^ 



15 AT—>oc R | 9 V-^oo ^5 



* to . (47) 

Si 



3 N >oo R 2 



r-l 



Knowing that lim^r »oo ^ is a finite non vanishing value, that 

limjv >ooN~ l Sn is a finite value, that J x 2 G^{x)dx = B N , and that 

K ' N = \ j ' ^ ' 4 U> N^ x ) dx ~ (J I x I 2 = j-fi Kn 



we conclude 



V^oo W 15 N-^oo N B l 9 iV^oo N N 



- 2 - lim ^^^0. 

3 iV^oc ^ ^3 

This completes the proof of our Lemma □ 
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Remark 3. IfV^ is a standard, confining, short-range and stable potential, at 
large N the entropy function Sn(v) = jr log Q (Nv, N) is an intensive quantity, 
that is 

Sw(v) ^ S N (v) . 

This is the obvious consequence of the well known fact that 

NS N (A d ,v)=N 1 S Nl (A d 1 ,v)+N 2 S N2 (Aiv) + 0(\ogN) (48) 

which is proved in textbooks \l^l] and which has also the important consequence 
summarized in the following remark. 

Remark 4. A consequence of equation fl^&| ) is that 

Q(Nv, Ni + N 2 , A( U Ag) = n(JViv, N x , A?) Q(N 2 v, N 2 , A d 2 ) 9{N) , (49) 

where 6(N) is such that [^(iV)] 1 ^ = 0{N l / N ) -> 1 for N -> oo. For two 
identical subsystems the potential energy is equally shared among them, with 
vanishing relative fluctuations in the N — ► oo limit. 

Remark 5. In the hypotheses of Theorem 1, V contains only short range in- 
teractions and its functional form does not change with N , i. e. the functions 
^ and $ in Definitions 3 and 4 do not depend on N. In other words, we are 
tackling physically homogeneous systems, which, at any N , can be considered 
as the union of smaller and identical subsystems. At large N, if a system is 
partitioned in a number k of sufficiently large subsystems, then the general- 
ization to k components of the factorization of configuration space given in 
Remark^holds. Therefore, the averages of functions of interacting variables, 
belonging to a given block, do not depend neither on the subsystems where they 
are computed (the potential functions are the same on each block after suitable 
relabeling of the variables), nor on the total number N of degrees of freedom. 

Lemma 8. Let {£i};=i,...,7v o,nd {yi}i=x,...,N be two independent sets of mutually 
independent non negative random quantities. Define X = YliLi x i an< ^ Y = 
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Yli=iVi- Let Y > for any realisation of the random variables {yi}i=i,...,N- Let 
(X) , (Y) denote the averages over an arbitrarily large number of realisations 
of the sets of random variables {xi}i=i,...,jv and {yi}i=i,... t N, respectively. 
In the limit N — > oo, it is 

W (Y) • 



Proof. After the Khinchin Theorem recalled below Definition El in the large 
N limit both X and Y are gaussian distributed random variables. Setting 
5X = X - (X) and 6(1/ Y) = 1/Y - (1/Y) we have 



Y / \Y I \ \Y 

Moreover 

6X S [ 1 ] \ < / 5Z 6 ( 1 



Y ) / ~ \ \Z / 
where Z = X if ((5X) 2 ) > ([5(1/F)] 2 ) or Z = Y if ((<5F) 2 ) > ((5X) 2 ), and 

w «(|)) = 1 - 2 < z >(^> + < z > , (f)- (51) 

Now, for a gaussian random variable Z such that (Z) > 0, we have 
1 \ 1 / 1 \ 1 



1+ <(^-<^» 2 > 



Z/ (Z) \1 + (Z-(Z))/(Z}/ (Z) L 3(Z)^ 

where all the terms with odd powers in the series expansion of 1/(1 + 6Z/ (Z)) 
vanish, and the even powers terms are powers of the quadratic term which is 
0(1/N), thus in the limit N -> oo 

^ 1 



Z/ (Z) 
Using Eq.(jn2} in Eq. (jBT|) we get 

5X5 (k))-- l + W) =0{l,Nh 

which, used in Eq. (j50|) together with Eq. (|52|) . leads to the final result. □ 
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5.2. Part B 

This part is devoted to the proof of the existence of uniform upper bounds 
as affirmed in the Lemma EJ 

We shall prove that the supremum on N and on v G 1$ exists of up to the 
fourth derivative of Sn(v). The proof of the existence of sup^ will be given by 
showing that the functions considered have a finite value in the N — > oo limit 
for any v G I~. The existence of the supremum on v is then a consequence of 
compactness j2y] of the set 1$. 

Remark 6. In what follows, the detailed proof is given for lattice potentials V^, 
however, in the fluid case the only difference is that the number of particles, 
interacting with a given one, is not preassigned. For this reason, in the fluid 
case, the number of particles within the interaction range of any other particle 
has to be replaced by its average. After the end of Section \5.2.°A more comments 
are given on this point. 

5.2.1. Proof ofswp N)ijeh \S N (v)\ < oo 

This directly comes from the intensive character of Sjy. □ 



where Q'(v, N) stands for the derivative of fl(v, N) with respect to the 
potential energy value v = Nv. 

The assumptions of our Main Theorem allow the use of the Federer- 
Laurence theorem enunciated in Section 0] and of the derivation formula given 
therein, thus 



5.2.2. Proof of mp N ^h ff(^) < °o 



By definition of SV we have 




1 Q'(v,N) dv _ n'(v,N) 
N Q(v,N) ' d/d ~ Q(v, N) 
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whence 

where ( stands for the configurational microcanonical average performed 
on the equipotential hypersurface of level v. 

Let us proceed to show that this derivative is bounded by a term which is 
independent of N. 

To ease notations we define 



X = 



w 



so that Eq. ()54j) now reads 



9Sn 



1 \ MC 

-m) 

A / N,v 



It is 



-Ax) 

X 



AV Jl^x&VdlV&iV 



and hence 



~Ax) 

X 



< 



|W|| 2 IIWI 



AV\ , JS =1 W(|W 
+ 2 ■ 



(55) 



(56) 



(57) 



|W|| 2 || W || 4 

where djV = dV/dq\ q % being the i-th coordinate of configuration space 

In the absence of critical points of V it is ||VV|| 2 > C > 0, thus we can 
apply Lemma |H1 where Y > is required, to find 



OS 



A' 



dv 



v) 



< 



< 



~A(X)) 
X I 



| AV 

\vv\\ 



/ N,v 



< 



Ax) 



/AC 



N,v 



+ 2 



N,v 



(I *v \)% v 



|w|| 2 >" c 



Nj 



+ 0[ N 



W|| 4 



+ 04, 



(\\w\\r N : v 

Consider now the term (| AV \)^ v - As the potential V is assumed smooth 
and bounded below, one has 



1 

N 2 



(I AV \)% v 



N 



i=l 



flC 



N 



< £<l 3^ l>R ^ N .m**(\%V\) 

' » i=l, ...N 



N,v 



i=l 
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As a consequence of Remark at large N (when the fluctuations of the aver- 
ages are vanishingly small) maxj = i ) .. ) jv(] d^V \)n v does not depend on N. The 
same holds for (I &V&V& i V iV? and max 1=1 N (I &V&V&V \)Z . 
We set mi = max i=lj .. i7V (| d\V |)^ andm 2 = max !J=lr . iW (| ^FflP-VW |>^ 
Let us now consider the terms (|| Vl A || 2n )^ c t; for n = 1,2. One has 



N \ i 10 N 



IwilX = J2(9,vf = £ < ( a,ioy;„ > w mm <(av)'>;. 



i=l / n v i=l 




£ <(s^) 2 (%v) 2 >;; 

By setting c x = mm i=lj „ tN ((diV) 2 )^^ and c 2 = min ij= i v . iA r {{diVf (9jV) 2 ) Nv 
we can finally write 



a /at 



where n p is the number of nearest neighbors. It is evident that in the limit 
N — > oo the r.h.s. of the equation above tends to the finite constant m\jc\. 
The upper bound thus obtained ensures that 8up_^ ee j s |%jf (^)| < oo. □ 

Remark 7. Notice that, in the fluid case, the computation of quantities like 
{{diV) 2 )^ v or (\df i V\) > ^ v involves an a-priori unknown number of neighbors 
of the i-th particle (we say that a particle is a neighbor of another one if the 
distance between the two particles is smaller than the interaction range of the 
potential). However, the requirement that V is repulsive at short distance, 
so that clusters of an arbitrary number of particles are forbidden, guarantees 
that each particle has a finite average number of neighbors. Thus, averaging 
quantities like the above mentioned ones yields N -independent values. 

In order to extend to the fluid case the proofs of uniform boundedness of the 
derivatives of the entropy ( qiven throughout the present Section one has 
to interpret n p as the average number of neighbors of a given particle. 
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Remark 8. Notice that the above computations show that 



lim \ - 



const < oo 



N,v 



which follows from the boundedness of \{A(x)/x}\- 



5. 2. 3. Proof of sup N>veh 



d 2 S N / — \ 



< oo 



The second derivative of Sjy can be rewritten in the form 



d 2 S 



N 



dv 2 



v) = N 



£l"(v,N) f£V(v,N) 



tl(v,N) \£l(v,N) 



(59) 



or, by using the same notations as before, 



d 2 S 



N 



~A(x) 

X I N,v 



IXC 



(60) 



again we are going to show that an upper bound, independent of N, exists also 
for this derivative. In order to make notations compact, we define 

V 



||W| 

N 



for any h ± ,h 2 , ^(h) . ip(h 2 ) = ^ 4>i(hi)4>i(h 2 ) 

i=i 

whence simple algebra yields 

±(V)-±(x) = X 2 M l -x^V , 

f{V) = ± {±{V)) = -±{V) ■ l(x) + X 2 AV 

TV 

mAV)) = x 2 div- x 2 MV)Y,MV)d 2 k v 



k=i 



N 



Mx) = -x 3 J2 d i v ^(v) 



3=1 



N 



k=l 



MUV)) = x 2 div- x 2 ^(v)J2MV)dlv 

1>i (d 2 r V) = xdtjrV 
i d l V ) = xdf j3 V 



(61) 
(62) 

(63) 
(64) 

(65) 

(66) 
(67) 
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where M l = V(W/||W||) 
notations we have 



-N ■ (mean curvature of £„). With these 



A 2 ( X ) = A (A( X )) = A U(V) ■ l(x) + X 3 AV) 



x - 



and thus Eq. (J5Uj) now reads 



d 2 S 



dv 2 



< N 



+ N 



A(X) 
X 





~/A(x) 


N,v 


\ X 


A(X)\ 






X J 


/ N,v 





1 2 



By using the relations (|61j) - (j67|) . the term -A (x) is rewritten as 



A(X) 



X 



-V> (-VWx) = -VW • i>{x) + x 2 av 

2 X M X - x 2 AV 



||W|| 2 l|V\/|| 4 
Now we consider the following inequalities 



iwi 



< 



< 



< 



||w|| 4 

M ; ^=i<i^i/^i>;; 



iV n p m 2 
c 2 N 2 



|w||%; 



(68) 



(69) 



(70) 



+ ° [ N~ 2 



(71) 



where n p is the number of nearest neighbours, and again 
m 2 = maiy^ (| ^V^V^V |>^. 

As m,2 keeps a finite value for liniAr^oo, the l.h.s. of equation (J7T|) vanishes 
in the iV — > oo limit. 

Thus, the larger iV the better the term (%) is approximated by £ = 
Eii<^/||W|| 2 = where £ = flgV/||VV|| 2 . Here we resort to the 

LemmaEland replace the microcanonical averages by "time" averages obtained 
along an ergodic stochastic process. Each term £j, for any i, can be then 
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considered as a stochastic process on the manifold E„ with a probability density 
Ui(£i). In presence of short range potentials, as prescribed in the hypotheses 
of our Main Theorem, and at large N, these processes are independent. 

By simply writing £ = Y^i=id = l/N^2 i=1 N£i, we are allowed to apply 
Lemma which tells us that the the second moment B' N of the distribution of 
£ is such that limAr^oo N B' N = c < oo. 

The first term of the r.h.s. of (|69|) is the second moment of - A (x) multiplied 
by N, this term, in the light of what we have just seen, remains finite in the 
N — > oo limit. 

Then we consider the second term of the r.h.s. of equation (|69p. This can 
be computed with simple algebra through the relations f)61ti67|) to give 



l(V)-± = 8 X 4 {(±(V)-±(V))) 2 - A X \±{V)\±{V)) 




N 



- 2 X \±(v)-±(v))/\v + x *Y.^ { y) d h v 



»J=1 



N 




(72) 



i,j,k=l 



where 



(i>{V);^{V)) 



(73) 



Wll 2 



(74) 



||W|| 2 

||VV|| 



(75) 



(76) 



Wll 3 
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The same kind of computation developed for equations (fTTj) gives 

nc N^nhn^ / 1 \ 



i a \ llc N 2 nlm 5 / 1 \ 

N { x \ t iV)\ t {V)))l v < -^ + 0^-) (78) 

N 3 n p m 6 „ / 1 



iV< X 4 (i(^);i(V))AV>- < ^-2ip + Ol^) (79) 



AT \ ^ c 



nU^W) ^^ + °(^) (80) 

I N \ MC N 2 n 2 m / 1 \ 

# ( x 3 E } < -— ^ +o(—) (si) 

\ i,i,fc=i / N)V 3 v y 

where, resorting again to the argument of Remark EJ we have defined the 
following quantities independent of N 

m 4 = . {{d i VdlVd j V){d k Vdl l Vd l V))^ v 

m 6 = . m^ N ((d i VdlVd j V)(dl k V))^ v 
max (diVdf^V)^ 



m 8 = . max ((d i Vd j Vd k V)d? jk V) tl ; v 



and 



c 3 =. mm ((d il Vf(d i2 V) 2 ---(d ie V) 2 ) 

ll,...,lQ=\,N 



N.v 



so that the r.h.s. of Eqs. (fTTJj) and (JHUj) have finite limits for N — > oo, while 
the r.h.s. of (J77|), (J7EJ) and (jHIJ) vanish in the limit N -> oo. 

In conclusion, since the ensemble of terms entering equation (JoUj) is bounded 



above, we have sup^g^ 



< oo. □ 



dv 2 

Remark 9. Notice that the above computations show that 
lim N (ip(V) ■ ip (^X^) \ = const < oo 



N^cc \~ - \ X J 
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5.2.1 Proof of sup N ^ h ^(v 



< oo 



The third derivative of Sn can be expressed as 
d 3 S N 



dv 3 
N 



(v) 



2 j n"*(v,N) 



z Q"(v : N)Q'(v,N) + 2 



Q(v,N) (tt(v,N)) 2 
or, by using Federer's operator A, 



n\v,N) 
n(v,N) 



d 3 s 



dv 3 



V) 



^ 3 (x)^ c 



N 2, / x _ 3 



/LtC 



^ 2 (x)y c M( X ) \^ j M( X ) y 

N,v \ X / N,v \ X / N.v \\ X / N,v, 



where 



X 



X - - V X / 



Ax) A 



x \ x J - - v x 



x x _ _ 



|W|| 2 

By substituting the expressions ([83 J) -([55 ]l into the r.h.s. of equation 
get 



d 3 S 



dv 3 



V] 



< N 



+ 3N' 



+ N 2 



ip(V)-ip[ ip(V)-ip 



x 



x 



N,v 

A X 



(83) 
(84) 
(85) 
, we 



v. 



/L1C 



AT,tA \ A 



X J \\ X J I .y.r 



3\ 



(86) 



By explicitly expanding the first term of the r.h.s. of (J86j) more than 30 terms 
are found. Nevertheless, these terms are similar or equal to those already 
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encountered above and, consequently, their iV-dependence can be similarly 
dominated as in the inequalities (j77H81j) . 

Consider now the second term of the r.h.s. of equation (jHEJ)- If we put 

using equations (j^ZJ) and (J72J) we can write 

N N 

t=l j=l 

Then 

(^^■*m>;-(^»-*(^)>; 

= E - W^<Pi>W • (87) 

Let us consider the terms, in the last sum, for which i and j label sites which 
are not nearest- neighbours j^. The corresponding expressions of a» and pj 
have no common coordinate variables. Thus, when computing microcanonical 
averages through "time" averages along the random Markov chains of Lemma 
IH1 we take advantage of the complete decorrelation of dj and pj so that 

for any i,j s.t. < i, j < N, then (a&j)^ - (a^^Jpj)^ = 

(where )i,j{ stands for i,j non nearest neighbours) which simplifies equation 

flEZD to 

(hi) 

< N n p max ((a iPj )^ - {a^ {p^ ) . 

(m) 

Now, equations (J58|) and (j77H81J) imply 

for any i,j s.t. < i, j < N, lim N 3 (a i p j ) l ^ v < oo 

N — >oo 

while equations (joTj) and (J72J) imply 

for any i,j s.t. < i, j < N, lim N 3 (a^ (pj)^ < oo , 

N — too 
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where stands for i,j nearest neighbours. Thus, the second term in the 
r.h.s. of equation (J86|) is bounded independently of N in the limit iV — > oo. 
The third term of the r.h.s. of equation (jHEj) is smaller than the third moment 
of the stochastic variable A(x)/x (multiplied by A^ 2 ). As we have already 
seen, we can rewrite A( X )/x = (U N ) Ei=i N9f i V/\\VV\\ 2 to which LemmaEI 
applies thus ensuring that the third moment C' N of the distribution of A(x)/x 

is such that \im N >oc A^ 2 C' N = 0. 

Finally we are left with a finite upper bound of the l.h.s. of equation (}86|) 
in the N oo limit. □ 

Remark 10. Notice that the computations above show that 

Jim N 2 (iAV) ■ V> (±{V) ■ ± ) ^ = const < 00 • 



5. 2. 5. Proof of swp N ^ h 



d i S N (x* 



< oo 



The fourth derivative of Sn{v) is given by the expression 



9 5jv /_>. 



Q iv (v,N) n"'(v,N) Q'(v,N) 



n(v,N) 

+ N 6 <J 12 



(n(v,N)y 



n"(v,N) 
n(v,N) 



, f , n n"(v,N) (tt(v,N))' 



6 



Qf(v,N) 



(n(v,N)y " \ to(y,N) 

Again we make use of the Federer operator A to rewrite it as 
'A\x)^ C 



9 Sn 
~9^ {v) 



N 



- N s 




QN 



X I N,v 

X I N,v 

., Mx) w x ' 



flC 



1( ^!My c /a(x)v 



fMC 



r2( my /my 

X / N,v \\ X I N,v, 



X 



X 



N,v. 
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where, after trivial algebra, 

m = (m) 1 +6 (m)\ (v) .jm) 

X \ X J \ X J V X J 

+ 3 ^(v) . ip f~^ry J +4^^(v) .i>(i>{v) . v ^ x j 



To make the notations more compact we use 



X 



v = tj)(y) . i> 



Ax) 
x 



w = il>{v) . i) ( . v ( 



so that, using again equations 



, we obtain 



d 4 S 



dv A 



N f -s 
V) 



< N 3 

+ 3iV 3 
+ AN 3 

+ QN 3 
+ N 3 



4\ MC 



39) 



jltC 



-4 - 



2 \ 



A'.i 



Consider the first term of equation (|89|). It is an iterative term already 
considered for the third derivative. This term stems from the application of 
the operator ip(V) ■ ip(-) to the term W which in its turn stems from the 
application of the same operator to the term V. The effect of this operator 
is to lower the N dependence of the function upon which it is applied by a 
factor N (what is simply due to the factor 1/|| W|| 2 ). Deriving with respect 
to v brings about a factor N in comparison to the derivation with respect to 
v, therefore the first term of equation is of the same order of N 2 (W)^„ 
and consequently, according to the Remark [TUJ it has a finite upper bound 
independent of N in the limit N — ► oo. 
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Consider now the second term of the r.h.s. of equation (|89[) . The Remark 
01 ensures that lini/v >oc N (V)^ v < oo. Moreover, after Lemma 

2 



(90) 



Consider now the third term of the r.h.s. of equation (|89j) . The Remarks El 
andfTUlent 
Lemma [7| 



and [TU1 entail limjy >oo(A)'^ v < oo and limjy >00 N 2 (W)^„ < oo. Thus, after 



lim^iVl (<W - (W)K„>^) < oo , 



whence 



lim N 3 

N — >oo 



A - (A)n,v)n, v 



lim N 3 

N — >oo 



< oo . 



(91) 



Consider now the fourth term of the r.h.s. of equation (|89p. If we write 



1 N 1 N 

i=i i=i 



with Oj and pj terms of order 1, we have 



N 3 



ftc 

'N,i 



V S ((^-Wi?*) ( a i-( a i)^«) {pk-{Pk)^ v 



i,j,k=l 



N,v 
lie 
N,v 
/ic 



where )i, j, fc( means that at least two of the three indexes refer to non near- 
est neighbours sites, whereas k) means that the three indexes are nearest 



neighbours. If k are such that k( then at least two of the three terms 
a,i, cij and pk have no common configurational variables. The microcanonical 
averages are again estimated according to Lemma El through a stochastic pro- 
cess on the configurational coordinates. The random processes associated with 
a*, dj and pk are thus completely decorrelated and one has 

for any i,j,k, s.t. )i,j,k(, 

( (ai - (a*)^) (a 3 - (p k - (pa>^„) )^ = . 

Now, if we consider i, j, k such that k), the three terms cij, aj and pk are 
certainly correlated but we notice that there are only Nn 2 terms of this kind. 
Thus we have 

< max | (ai - (a^) , (p* - (Pfc)^„) } • 

Since the terms Oj and pk are of order 1, the largest term of the preceding 
equation is independent of N, we have thus found the upper bound of the 
fourth term of the r.h.s. of equation (JHHJ). 

Finally, the last term of the r.h.s. of equation (J8§j) is the fourth cumulant 
of the stochastic variable A(x)/x (multiplied by N 3 ). As already seen above, 
we write A(x)/x = V^Z^i -^^V/ll VV|| 2 so that Lemma [3 applies and 
ensures that the distribution of A(x)/x nas a fourth cumulant K' N such that 
lim^oo N 3 K' N = 0. 

The ensemble of the upper bounds thus obtained yields the final desired 
result. □ 

6. FINAL REMARKS 

To conclude this first paper, some comments are in order. 

Remark 11 (Domain of physical applications). Notice that the require- 
ment of standard, stable, confining and short-range potentials applies to a 



broad class of physically relevant models. In fact, the interatomic and inter- 
molecular interaction potentials (like Lennard-Jones, Morse, van der Waals 
potentials) which are typically encountered in condensed matter theory, as well 
as classical spin potentials, fulfil these requirements. 

Remark 12 (Sufficiency conditions). Notice that the converse of our Main 
Theorem is not true, in other words there is not a one-to-one correspon- 
dence between any topology change of the energy level sets and phase tran- 
sitions. In fact, there are systems, like the Fermi- Pasta-Ulam model described 
by Vw(g) = YuiLi _ <??) 2 + f (Qz+i —QiY which, for fixed end points, has no 

critical points and no phase transitions, whereas, for example, a one dimen- 
sional lattice of classical spins (or of coupled rotators) described by the potential 
function Vj^(q) = YliLii^ ~ cos (.Qi+i ~ Qi)] has many critical points \l A] so that 
both families {H v } v( z^ and {M v } v& ^ undergo many topology changes, but, since 
no phase transition is associated with this potential, none of these topology 
changes corresponds to a phase transition. Note that this is not a counter ex- 
ample of our Main Theorem (which would require to find a system undergoing 
a phase transition in the absence of topology changes and within the domain 
of validity of the Theorem), it just tells us that the loss of diffeomorphicity of 
the {E v } ve ^ and, equivalently , of the {M„}«£i at some v c , is a necessary but 
not sufficient condition for the occurrence of a phase transition. 

Remark 13 (Relevance of topology changes for phase transitions). 

In order to prove that our Theorem is relevant to statistical mechanics, and 
in particular in order to really link the phenomenon of phase transitions to 
a topology change of the configuration space submanifolds M v , in paper II we 
work out an analytic relation between configurational entropy S(v) and the 
Morse indexes of the submanifolds M v . Such a relation is formulated within 
another Theorem (enunciated also in the Introduction of the present paper) 
which unveils why the differentiability class of S(v), in the N — > oo limit, can 
be lowered from C°° to C 2 or to C 1 only by a suitable energy change of the 
Morse indexes (hence of topology change). Loosely speaking, in the context of 



45 



our topological approach, the Theorem proved in paper II plays an analogous 
role to that played by the Lee- Yang circle Theorem ]2R] within the context of 
the Yang-Lee theory of phase transitions. 
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